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^ Phase insensitive optical amplification of an unknown quantum state is known to be a fundamentally 

Ch noisy operation that inevitably adds noise to the amplified state [1-5]. However, this fundamental noise 

|JZ penalty in amplification can be circumvented by resorting to a probabilistic scheme as recently proposed 
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and demonstrated in refs ||6H8l. These amplifiers are based on highly non-classical resources in a complex 

interferometer. Here we demonstrate a probabilistic quantum amplifier beating the fundamental quantum 

(— I limit utilizing a thermal noise source and a photon number subtraction scheme [9|. The experiment shows, 

CL| surprisingly, that the addition of incoherent noise leads to a noiselessly amplified output state with a phase 

-<-j uncertainty below the uncertainty of the state prior to amplification. This amplifier might become a valuable 

S quantum tool in future quantum metrological schemes and quantum communication protocols. 

^ Besides being the subject of a fundamental dicussion going back to Dirac [ 10|, the measurement of phase 

, ^ is at the heart of many quantum metrological and quantum informational applications such as gravitational 

wave detection, global positioning, clock syncronization, quantum computing and quantum key distribution. 

S^ In many of these applications the phase is most often imprinted onto a coherent state of light and subse- 

\Q quently estimated using an interferometric measurement scheme. Such a phase estimation process [11] is 

O however hampered by the fundamental quantum noise of the coherent state which plays an increasingly 

devastating role as the excitation of the coherent state becomes smaller. Small coherent state excitations 

and associated large phase uncertainties are typical in real systems such as long distance coherent state 

Q communication and lossy interferometry. 

O To reduce the phase uncertainty and thus concentrate the phase information, the state must be amplified 

^ noiselessly. This can be done probabilistically using either a highly complicated interferometric setup of 

J> single photon sources |^-'8l, a sophisticated sequence of photon addition and subtraction schemes ll9l[T2l 

KA or a very strong cross-Kerr nonlinearity [13]. However, as we show in this letter, it is possible to amplify 

H the phase information noiselessly without the use of any non-classical resources or any strong parametric 

interactions. Remarkably, the supply of energy in our amplifier is simply a thermal light source. 

A schematic diagram of the probabilistic amplifier [9 ] is shown in Fig. [T^. It is solely based on phase 
insensitive noise addition and photon subtraction. To explain in simple terms why the addition of noise can 
help amplifying a coherent state, we consider the phase space pictures in Fig. [TJ5. The addition of thermal 
noise induces random displacements to the coherent state, thus resulting in a Gaussian mixture of coherent 
states; some with excitations that are larger than the original excitation and some with smaller excitations. 
In the photon subtraction process, the coherent states with large excitations are probabilistically heralded, 
thereby rendering the state in a mixture consisting of the most excited coherent states from the original 
Gaussian mixture. As illustrated in Fig. [l]?, the resulting state is amplified and possess a reduced phase 
uncertainty. 

The probabilistic photon subtraction procedure can be approximated by a largely asymmetric beam 
splitter combined with a photon number resolving detector, PNRD (see Fig. [T^). A small portion of the 
displaced thermal state is directed to the photon counter and when a pre-specified number of photons is 



detected the transmitted state is heralded. Such an approach for photon number subtraction has also been 
employed for the generation of coherent state superpositions lfT4l [T5]| . However, in contrast to previous 
implementations that were limited to the demonstration of two-photon subtraction 1.16.1 . here we subtract up 
to four photons. 

To elucidate the function of the amplifier, theoretically, we consider the amplification of a small ampli- 
tude (\a\ <^ 1) coherent state which can be approximately described in the two-dimensional Fock space: 
\a) f« |0) + a|l). Since the amplitude is small, the canonical phase variance of this state is to a very good 
approximation given by IITtII 

Vc-t\ (1) 

This variance represents the fundamental uncertainty in estimating the phase of the coherent state when a 
hypothetically ideal phase measurement is employed lITSl . The aim is to produce an amplified state with a 
phase variance reduced with respect to the coherent state variance in ([T]l, thereby concentrating the phase 
information. If a conventional phase insensitive amplifier is used to amplify the coherent state, the resulting 
variance is larger than ([T]l (see Apendix B). On the other hand, if our amplifier is employed with weak 
Gaussian noise addition followed by single photon subtraction, the resulting state is ||9| 

P ^ I \2^M \a\ \2m [I«I'|0)(0| + ^th(lO) + 2a|l))((0| + (l|2a*)] (2) 

with the canonical phase variance 

^^'"'^^^ (' + €)-'■ <^' 

where it is assumed that the average number of incoherently added photons is A'th ^ 1- We quantify 
the performance of the amplifier by the normalized phase variance, F = V^^^ /Vc which is smaller than 
one for a noiseless operation. For the above approximative example, if |ap <^ Nth, we find that the 
normalized variance approaches T = 1/4. Another parameter that will be used to evaluate the amplifier 
is the gain g = \/3\/\a\ (where /3 is the average amplitude of the output state), being 51 = 2 for the above 
example. We therefore see that by simply adding a small amount of noise to the input state followed by 
single photon subtraction it is possible to create an output state with twice the amplitude and with a reduced 
phase variance. 

Based on a more general model (as presented in Apendix C), in Fig. |2] we plot the normalized phase 
variance and the gain as a function of the average number of added thermal photons, where M denotes the 
threshold for the number of photon subtractions. The figures illustrate three interesting aspects: the phase 
noise reducing operation works even when the parameters go beyond the simple approximation considered 
above, the effect of the amplification improves with subtracting more photons, and the amount of noise 
must be relatively large for the amplification to work well. These aspects also follow the intuitive picture 
discussed above: Large noise addition will partly displace the coherent state in a radial direction in phase 
space with a correspondingly large magnitude, and the states with the largest amplitude (associated with 
strongly amplified states) are heralded by high photon number subtractions. We also note that the amount 
of added noise that minimizes the canonical variance depends solely on the magnitude of the input coherent 
state but not on its phase. The amplifier is thus capable of concentrating an unknown phase of a coherent 
state. 

A laboratory implementation of the amplifier is depicted in Fig. [TJ;. Our source is a grating stabilized 
CW diode laser operating at 809 nm with a coherence time of 1 /is. The laser output is spatially cleaned in 
an optical fiber, and subsequently split to serve as a Local Oscillator (LO) for homodyne detection and as an 
auxiliary beam for state preparation. We use a pair of electro-optical modulators for the preparation of the 



displaced thermal state (corresponding to a coherent state with added thermal noise) in a polarization mode 
orthogonal to the polarization of the auxiliary mode (see Apendix A). The duration of the prepared pulses 
is 800 ns. A portion (20%) of the prepared state is tapped off in an asymmetric beam splitter and measured 
with an avalanche photo detector (APD) operating in an actively gated mode such that the dead time (50 ns) 
is much shorter than the pulse duration. This means that the APD can be used as a single photon counter 
provided that the mean number of photons in the detected pulse is very small. The transmitted part of the 
state is passed on to the homodyne detector where it interferes with a phase controlled LO. This provides 
quadrature measurements of the emerging states under any phase space angle. The measurement outcome 
is sent to the computer where it is postselected according to the result of the photon counting measurement. 
Based on the resulting data points we reconstruct the density matrix of the heralded state using a maximum 
likelihood algorithm 1 19, 20|. We correct the data for the inefficiency of the homodyne detector in order to 
reconstruct the actual input state and the amplified output state (see Apendix A). 

From the density matrices, we construct the Wigner functions for the input state and the amplified 
output states for different photon number subtractions as illustrated in Fig. [3] Here we consider an input 
state excited along the amplitude quadrature axis with {X) = 0.431 thus |ap = 0.186 and a thermal noise 
addition corresponding to (A^th) = 0.15. The amplification factor for this experiment is summarized in 
Fig. [4j;. We also reconstruct the phase distributions (see Apendix B) for different subtractions, the results 
of which are shown in Fig.|4^. We clearly see that as the number of subtractions increases, the distribution 
becomes narrower, and thus the phase information is concentrated. These results are summarised in Fig.|4j3. 

To optimise the performance of the amplifier - that is, to minimize the phase variance - the amount 
of added thermal noise should be chosen appropriately with respect to the input coherent state amplitude. 
Furthermore, we note that having detailed information about the input alphabet, the structure of the noisy 
displacements can be tailored accordingly, thereby drastically reducing the amount of energy used to drive 
the amplifier. For example if the input is a phase-covariant coherent state alphabet, the optimized structure 
of the noisy displacement is also phase covariant. Such tailoring of the displacements as well as applications 
of the amplifier will be interesting directions for future research. Finally, we note that the noise addition 
process can be also carried out with a linear amplifier. Such an approach will not only add thermal noise to 
the input state but will also displace it coherently in the prefened direction, thereby further concentrating 
the phase information. 

In summary, we have reduced the phase uncertainty of a coherent state of light through noiseless proba- 
bilistic amplification. In contrast to previous approaches to noisefree amplification, the amplifier is neither 
based on an ample supply of nonclassical resources nor on strong parametric interactions, but solely on 
Gaussian noise addition and photon counting. Due to its pivotal properties such as simplicity and robust- 
ness, we expect that this approach to probabilistic noisefree amplification will be of interest to a large 
variety of experiments and protocols involving phase estimation such as quantum metrology and quantum 
communication. 



APENDIX A 

The experimental setup in Fig. [T]: is described in the following. The laser source is monitored, to assure 
quantum-noise-limited signal states. The states are prepared by two electro optical modulators (EOM) and 
a half wave plate (HWP). The modulators displace the signal state (S) using the orthogonally polarized 
auxiliary oscillator mode (AO), which is relatively bright 1211 . After the modulators' calibration, we can 
displace the signal state to any coherent state with a maximum photon number corresponding to the mean 
photon number in the AO mode rimax = |aAOp- The signal mode can be chosen to be in an arbitrary mixed 
state, provided that the state's P-function is positive. We can therefore generate the displaced thermal state 
applying a suitable modulation sequence to both EOMs. The state is modeled with a finite set of (more than) 
10^ coherent states, randomly picked from a 2D normal distribution. The random modulation sequence is 
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FIG. 1 : a) Principal setup: Noise is incoherently added on an input state. Subsequently a part of the beam is tapped off 
and measured by a photon number resolving detector (PNRD). Results of that measurement that surpass a specified 
threshold M, herald the output state, b) Principal operation in phase space: A coherent state (yellow) serves as the 
input. The dashed line indicates the phase variance. Thermal noise is added to this input state resulting in a displaced 
thermal state (green). The output state (blue) is reshaped and the resulting phase variance is reduced compared to the 
input state, c) Experimental setup: An external cavity diode laser with fiber mode cleaning (FMC) acts as the source 
for the experiment and is split into a local oscillator (LO) and an auxiliary oscillator. The signal is prepared by a 
combination of two electro-optical modulators (EOM), a half-wave plate (HWP) and attenuation (Att.). A polarizing 
beam splitter removes the auxiliary oscillator. Part of the signal is tapped by a 80:20 beam splitter and coupled via a 
multimode fibre into the photo number resolving detector. This measurement is conditioning the output state, which 
is characterized by a homodyne measurement. The phase of the LO is controlled by a piezoelectric transducer (PZT). 



varied and repeated throughout the measurement. 

The tap beam is focused into the fibre coupled PNRD (for details see f^Sl. '231). The transmitted part 
of the beam is sent to a homodyne detector, which measures the signal with a continuously-scanned local 
oscillators' phase. The scanning frequency is chosen to be 21 niHz, leading to an effective phase drift of 
only 1.6 mrad within the modulation sequence. This value is negligible from an experimental point of view 
so that the LO's phase is considered constant within a single modulation sequence. The phase angle needed 
in the tomography was estimated with a series of phase calibration signals prepended to the modulation 
sequence. The main source of error in the setup is the drift of the modulators. To compensate for this drift, 
the calibration point was continuously adjusted. 

As the homodyne detector (HD) is not a part of the phase concentration scheme itself but only imple- 
mented in the setup to prove the effect of the scheme, we do not want to take its imperfections into account 
in the analysis. We assume therefore perfect detection. The amplitude of the coherent input state |a| is then 
inferred from the measured mean photon number values in the imperfect PNRD and the ideal HD 

1 , 
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|aPNRD| 



(4) 
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where the PNRD's quantum efficiency was calibrated to be r/PNRD = 0.63 ± 3% using the overall quantum 
efficiency of the HD. This procedure is preferable as it does not demand an accurate knowledge of the input 
coherent states amplitude, the splitting ratio and the losses in the HD. 



APENDIX B 



Each measurement procedure devised to estimate a phase of a given quantum state can be characterized 
by a real semi-definite matrix H. The actual probability distribution of the estimated phase value can 
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FIG. 2: Theoretical gain and normalized phase variance vs. the mean number of added thermal photons: a) The 
canonical variance normalized to the corresponding variance of the input coherent state, b) The gain generally grows 
with the number of photons in the added thermal noise and the threshold M of added thermal photons for an input 
coherent state of amplitude \a\ = 0.48. 

be obtained as P{6) = Tr[pF{e)], where F{e) = l/27r Xl^^^g exp(i6'(m - n))iJ^„|m)(n| 117]. To 
characterize the quality of phase encoding by a single value, we can use the Holevo phase variance Il24l 
V = |/u|~^ — 1, where /i = {exp{i6)). The fundamental limit of the phase estimation is obtained for the 
canonical measurement, in which the operator F{9) projects onto the idealized phase state Yl'^=o e*^"i?^) 
and H = 1. The canonical phase variance is therefore given by 
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FIG. 3: Reconstructed Wigner functions: Wigner functions that were reconstructed from experimental data for the 
input state and the heralded state for different thresholds M. For each experimental reconstruction, the mean values 
and standard deviations for the X and P quadratures are given with the corresponding measured success probability 
and canonical phase variance. 
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FIG. 4: a) Phase probability distribution function derived from the experimental data (solid lines) for the measured 
coherent, the thermal and the conditioned states. Corresponding theoretical functions (dashed lines) were calculated 
for states fitting to experimentally derived parameters, b) The canonical phase variance deduced from the experimental 
data (light blue) and corresponding theoretical values (dark blue) calculated for states fitting to the experimentally 
derived parameters. The inferred input coherent state serves as the reference value, c) Gain for the input coherent 
state for different thresholds M. The error bars in b) represent the statistical deviations over many different realizations 
of the experiment. 



APENDIX C 

In the following we theoretically describe the amplification process and explicitly calculate values of ^ 
for various states. For a canonical measurement of an initial coherent state with amplitude a, we find |[9]| 
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It can be seen that the fundamental limit on phase estimation given by the canonical variance is fully given 
by the mean number of photons of the state. When the coherent state gets amplified by a parametric device 
with gain G, the result is a displaced thermal state state with a' = \fGa, A'th = 2G — 2, and |,9J 
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There is another way of obtaining fi, which will be useful later. At its core stands the Glauber-Sudarshan 
representation ESll . which allows any quantum state to be expressed as 



p= m3m{l3\d^l3, 



(8) 



where the $(/3) is its P-function. For a coherent state it is a delta function 5(/3 — a), while for a coherent 
state with thermal noise it is a Gaussian function 
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where Nth denotes the mean number of thermal photons. Since ^u is a linear functional of the density 
operator, we can, for any state whose P-function is known, obtain its value directly from Q and (|6]l. For a 
displaced thermal state with (|9]) we arrive at 
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Here we have defined the function In{A, B) of a positive integer n, a complex number A = Ar + iA-i and 
a real number B, as 

X„(A, B) = j{Pr + imPl + PlTGpMr, B)Gp^{Au B)dMPi 

where Qp {A,B) = exp(— (/3. — A)"^ / 2B) / ^/IttB is a Gaussian distribution of a variable /3. with mean 
value A and variance B, and A^fc(j4, B) stands for its A;-th moment. The moments can be easily evaluated 
with the help of the moment generating function 
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or a recurrence relation 



Mk+i{A, B) = AMk{A, B) + BkMk-i{A, B) 



(13) 



withA^o(^,-B) = 1- 

The process of amplification can be simplistically viewed as applying a number of annihilation opera- 
tors to the density operator of a sufficiently thermalized coherent state with pc given by ^ and (|9]). The 
amplified state is then proportional (up to normalization) to a^ pca)^ , where a and a) are annihilation 
and creation operators and /i can be easily found. However, although this approach gives a reliable ap- 
proximation for the idealized version of the tapping and measurement processes and provides a worthwhile 
insight into the limits of the method, the realistic procedure with a finite reflectivity of a beam splitter and 
inefficient photon detection is better served by the P-function treatment. In this picture, the amplified state 
can be expressed as 
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where $(/3) is the P-function of a the coherent state with added noise ([9]), Ps is the success probability and 
V{I3) represents the probability of a successful post-selection event for a coherent state |/3). For a realistic 
M photon measurement implemented by a beam splitter with transmissivity T and a photo-detector with 
quantum efficiency rj, the 'P(/3) can be expressed as 



Vnm = (V^(i-r)/3|n|Vr?(i-r)/3), 



(15) 



where we have modeled the detector inefficiency by a virtual beam splitter with transmissivity ry, and fl = 
1 — X]a;=o I ^) (^I represents the positive detection POVM element of the ideal detector. We can then obtain 
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with Hi = NthT + 1 and E2 = Nt^T + Nttrjil — T) + 1. For the description to be complete we need 
to evaluate the success probabihty Ps = f ^W/VT)'Pu{f3/VT)d?l3/T. In a complete analogy with the 
calculations used previously we can express the success probability as 



Ps = l-- 
where H3 = iVth??(l - T) + 1 and 
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is defined analogically to ( 1 1 1. 
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